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4* e r p c £ he . A §~&r ~to r SioncJ Spring S^) . 



32 






k t( + Kfc* + ^ 1 3 7T£ 1 P,’ "2 



iz ^3 

+ c "t* h 

h pj p7 



32. 

ttg 



t, d/ P? + is. J T>i + A J>/ 

J>? P,* 1 !>/ 



or 



TIG D, D 2 D ? 



32 P 3* + ^2 P 3 f + ^ ^ ) 

Cc ") steps bebcuve <xS ^fifi'es /SprmgS. 



CO 
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Q) (*) f “ + = ( 50 ° x + 



L o c/jc I Xo 

~ iioo 'x. _ 4 ooo 

(b) oX x= 9 mm : , 

B xcuct F 9 = ^00X9+2 (?) = 535-8 N 

Approx imcu’ts Fg — | j o 0 x 3 — 4 OOO — ^ ^OO N 

Error 9 7 3 5 / 

(c) oX X — II mm: 3 

Exact Fj, ~ 500 X11+ 2 00 = 8 162 N 

Afifirc.xrma.-te F,, = Itoo x| 1 - 4000 = 8100 N 
Error =• + o »75 ’36 • 






= Co-nsi^t --(El') ; WfferenbiaAton of 1 ) 2' ve * 

Jtjo V f + f t dv - o 

^ —CeO 




ckcun^e in volume Udhen 


ma.cs’ movey 


£g, s • (Ez") and (e 3 ^) give 


4= KA 




V 


Force due to fress-ure 


change = . 


Spring constant of air Spring = ^ 



_ Jf 






(l3o) 



Equivalent Spring constant* in Jf'ffent cUrec-tfonS are 
. / V 5 ^7 \ . _ / ** fc 9 N 

61 ”(*5*6 + k 5 *7 + *6 *7 ) ' eZ 

\ + << 2 / * v + ** / 



= 



+ ) v Xj -r 

If -the force F move? by x, spring /tocated of 0^- undergoes cu 
JliS placement of = X cos 0^ (derivation as ,n problem M7)* 

Equivalence of potential. energy giVeS l^eg X = f j 

= 2 (* e ; “S 2 ®0 



t=i 



1-23 



PRELIMINARY VERSION 12/14/2010 



© 201 1 Pearson Education, Inc., Upper Saddle River, NJ. All rights reserved. This publication is protected by Copyright and written permission 
should be obtained from the publisher prior to any prohibited reproduction, storage in a retrieval system, or transmission in any form or by any 
means, electronic, mechanical, photocopying, recording, or likewise. For information regarding permission(s), write to: Rights and Permissions 
Department, Pearson Education, Inc., Upper Saddle River, NJ 07458. 




- (t) - ^ ^ “ fc* * ( 2 t ) + f ( £ ) — o 



or 



F = (^ + 4 ^ 



3 -2 + 

If denotes {-he u •‘v<=dje^d: 07 v^£cuy^ 

of- the C. Tke dlfS-cFlOn of F ad j?Oi'rA~ Cj 

uie ko^ve 

( 2 ) 

CO ° , - ryv d C 2 9* VC 

■ + 2 — ^ -+• = ^- f - y (^^)+( 3 ' #< ) 

O) 



- 
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Spfi’n^ Cor»sto-r \ t of a- -f >e\ica\ /-SJpn'n-g is 



'k = 



GcT 
S w P 



CO 



Assuming sfe.cc r modulus of 

CO ^iVes , for X> - o ' 2 . rr\ J dO 

1c = (7?-3 xio 9 ) (o -oo 
S (l o') (o-2) 3 



sfeel qa G = 11*3 G fieu, 
o-oo 5 m cuv-^d! jsj ~ (O^ 

= 77-V^I^ N/m 
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(A) 2) a.v>d d : jcLmefoi- fcotk belicoul Springs 

°f cl, be\icoX ip ring is: 

W = 7T r> mf) Nf ro 

where Z~ - specifi’c wei^Kt of moJteriaX of spring. 

For qj A\e.e I spring with - -?Q> -S’ ‘U N/ rn 3 > 

we ht is ( for - io) •• 

W* = TT X> (Tret 2 -) N* ^ ( 10 ) (76 '5 * 1 /) 

^ 4 

= l?-l2Sxi0^7r 2 DJ Z 0 2 ) 

For cly> o\x*-rr \ i r> w m /ijorirtj with 2.6-6 teN/xir?} "FAe 

Ivei^Kt ^ (" for nuwkr of furnr /V^, ) 

- 7T D QlJ* 2 ’ ) ^ 7T 2 D x //) 

= 6*65 * <0 TT p d 

E^Ucdln^ (z") curut (?)j 

1^*115 xio^ TT 2 ” D d *" = 6*6*3 *10 TV J) d 

or ^a> = l<? ‘ 1 2 * * = 28*75?^ 

6 *6 5* id 3 

^ Spring cans to- of a* ^elicoA ^Sprirv^ 1 S : 

ic = C ^ / (8 W D * ) 

For cb ^feel spring with & - ~l?3-3 G? co * 

\l a - (7?* 3 m c ) J § (' 10 ) D*} 

— 0-9.9 12. 5 x io ? P 3 

For oXu.tr> in um AprircJ with G — 26*2 G 

= C 2 6 ' 2 * {0 ) { S (ZS.7S?4) p* j- 

-- © • I / 3 7 x 1 0^ d ^ Z ' p ? 



ro 



co 



(5 ) 



rO 



*^ s ‘CO ( 6 ) (V\ d i oa 4 x tkaA ike 'Spring Coy^A^Jt of 

Steel sirring U °' 39 ‘ 2 ’</ 3 ? - S*“ 7 o ^6 times larger 



Tho^n t 


njf ; p, M spring 

ror^ i~ 1 iju/i im a n\z.Jur r-> cu ^ k i a r\ ia a /on h a 
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From Fro Hem U — f A u/Hh f- I* 4 £or air 

Let f> = 200 psl v z 

, (200)( |, 'O A 2- _A__ - o.2679 

*= 75 a/m = * V z 2 

Let di'cun eter of pi'stcno -A - 2 inch > A — "^T — 3>I4I£ 

= f< Z /o- 2679 = 36-S40S in 3 

Let |y = 2 . TncL ; 31 (z) V 3 ) - 4 -*? 4 2 9 inch 




1.36 



F=ax + bx 3 =2 (10 4 ) x + 4 (10 7 ) x 3 

dF * 

Around x : F(x) ~ F(x ) + I x* ( x — x ) 



When x* = 10“ 2 m, F(x*) = 2 (10 4 ) (10“ 2 ) + 4 (10 7 ) (10~ 6 ) = 240 N 



— | x -=a + 3bx 2 =2 (10 4 ) + 3 (4) (10 7 ) (10“ 4 ) = 32000 
dx 

Hence F(x) = 240 + 32000 (x - 0.01) = (32000 x - 80 ) N 

Since tbe linearized spring constant is given by F(x) = keq x, we have kgq = 32,000 
N/m. 
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F* = OJ I x- -f- } i- 1 , 2 . 



Sprimjs 1 in series : 

VJ - c > i S, fO 

W — Qsq 52 + y 2 . ^2. (2) 

W=^ e& S st (3) 

= 8, + S z (4) 



Solve Egf-O) cur»d (2) for £| cuiA 5*2, > 
r e s p eetiVely • Substitute the result 
in £^-C4) then in Eg,. ( 3> to 

£-i r>d ’ 
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/J'TX ic = —J - — > 8X10 N /m ' X - 6 : 

8 D 3 M 

KDbJf (TLi- 1 ) wkere f- w f>er unft volume 



- T i*JL = T 



SJ 2 ? 



w 



Z ^ 2 7T 2 b* NJ 2 J 5 



£ O • H- "2- 
j> - 76000 M/m 3 , 2= 9-Si m/sec 2 -, 



U Sin 3 G - 7 3 > I x to 9 n/ id 2 
— = 6, 8, 10 ; NJ = 10, 15. 20 ; J -0-4, 

Cl 



Vculues 



axe 



ccm jbuied . 




1.39 



•fc a_n el ^ j 

Co m b fna/fc ib n b , M = 10 txnj J = 2*0 im ? cl intf 

d 

° -\C - S - 4-6 0 & X io N/m CLt\cl ^ Hz , Ca-n be 

fczfe-e-n aJ> an ux.ce.pb coble Aesi'yci* 

Total elongation (strain) is same in each material: 



e 3 e a e 



where x is the total elongation. Equation (l) can be expressed as 



E= 



E, 



x 

_ _ 
E. x 



or cr. = 



v* = 



9 € 

E a x 



(1) 

( 2 ) 

( 3 ) 

( 4 ) 



Total axial force is: 



F = F 9 + F a = <x s A, + (T a A a 



( 5 ) 



where F s and F a denote the axial forces acting on steel and aluminum, respectively, and 
a and A, represent the cross-sectional areas of the two materials. Equating F to x 
where keq denotes the equivalent spring constant of the bimetallic bar, we obtain from 

Eqs. (3) to (5): 



F — kgq X 



E s x 



A* + 



E a x 



A a 



Eg Ag ( E a A a 

or k^ = — -j , — H ~p 



( 6 ) 
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o 








Let tee fen^te o f tkc 

SjfTirxj f>e i ■ S p r{r \? is 

unctfi ^orm eet ok Q-O- 
U/ken tfce enJl A °f 
Ike Spring i'S AiS 
by cun oumount yc cJ% 

^kou/h in tee figure.; 
tke /Spring is stretched 
by "the ctmouni X -A so tfkcufc Tke ■force. 

in tke (F^') is given ky 

r * - k) 

l“ke Component: ojl Ike fSprtr\j -force Rg oJI&y-lJ tee 
Jireclron of oc <5 given tg 



F* - F, e - F 



v* 



£ 



9 2 

z + X 



_ 



I r 
+ * 



-o 



Jl 



z -f- x z 



= k / 1 - 



u 



\fh 



X 



CO 



a ^ 

•f < 



Eg^cJtlav* (2) About fAcJ- 'tec Torcc - 

(in tee cUre.ctzcm') is nonlinear • x j- 
|ke reAa-&« is linecur, coedcA tVriFe 

CO 



Foe. = 



k x 



A C srrr^O^Ji^rryn or^ E^S • Q 2 - ) o-nru?( ti^edt 

ffc e ,^TXng * IS a- tu ^ 

Jte^en-X* cm tee x * 
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From -H\e 31‘ver* cta,-fco/, the force - fa rma/fct'on 

relation <?£ e Spring ir< obha-urULeA \ 



Tensile Torce ( F ) ^ 0 


too 2 50 3?o n$c> 570 






13 33 44 64 76 




O 


(cJrdSJy^e- in 







-j-j'xg foroe - i t* € I i o n *5 plo'FFfin in tfce 

figure below- Tke region ca -n be /Jeen to be 

« eo .r( 3 linear t*e ifri^ conS+o-nl" 3<'v«n by 



He = 






— 76 



5_Zi = 7.5 lV/-m trt =7 5-0 0 Nj/r 



F 
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d0 = 



T dx 



T dx 



GJ (80 (10 9 )) 1.5708 (0.1 - 0.05 x) 4 
The total angle of twist can be determined by integrating Eq. (1) from x=0 to 1 as: 

✓ \ i 



\ T dx 

0 = J 



J (12.5664 (10 10 )) (0.1 - 0.05 x) 4 



12.5664 (10 10 ) 

V ' 



dx 



But / 



dx 



o (0.1 — 0.05 x) 4 
dy 



( 1 ) 



( 2 ) 



1 v -0.05 J 

__ = L_ f (- °- 05 te) = - 20 f — rj 

\* 0.05 { (0.1 - 0.05 x) 4 0 (0.1 + y) 



J o (0.1 - 0.05 x) 4 °- U5 o (0. 

= 4.6667 (10 4 ) where y = — 0.05 x 

a T (4.6667) (10 4 j T ( 0 .3714 (10~ 8 )) rad 
Hence 0 - ^ l 

This gives k, = 2-6925 (10°) N-m/rad 

© The steel end aluminum hollow shafts can be treated as two torsional springs in parallel. 
For a hollow shaft, 
kr = ^ G - (D 4 - d 4 ) 

For the steel shaft, G = 80 (10 9 ) Pa, t = 5 m, D = 0.25 m, d = 0.15 m, and hence 
k _ 7T (8 (lo 10 )) ^ 0 .25 4 — 0.15 4 ) = 5.34072 (10 8 ) N— m/rad 
fa) For the ataSmm shaft, G = 26 (10 s ) Pa, t = 5 m, D = 0.15 m, d = 0.1 m, and 
hence 

k = 7T (26 (10 9 )1 (0.15 4 - 0.10 4 ) = 0.207395 (10 6 ) N-m/rad 

= kti + = 5.34072 (10 8 ) + 0.20739 (10 8 ) = 5.54811 (10 8 ) N-m/rad 

(£> ^ Wtik G{ = 26 (/0 9 ) Feu j jl ~ 5 7> - c -15 m and A - ° ’ 0 ** ^ , 
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k t - ^ ( 2 6 * (Vi* 4 - 0 . 0 5 4 ) = 0-255255 x \0 N 
2 32 (50 V ' 



U e y r +><tz = 5*3 lio-l 



r s- & 

2 xi o + 0-2552 55 xic ~ 5-595T915 xio 

N-m/ retd 



G d 4 



[1.44) For helical spring: k = 64 nR 3 



Spring l: k. - ^(ioxP t ° 1 ’ 388 ' 89 lb/b 



Spring 2: = 50.00 Ib/m 



(a) Spring 2 inside spring 1 (parallel): keq = k x + k 2 - 1,438.89 lb/m 

(b) Spring 2 on top of spring 1 (series): 



1 1 ^ 1 _ k 2 + kx 

“ kT k 2 



ki k 2 



which gives k^ = 48.2625 lb/in. 



L45) For oj be! iced spring V = ^ — 










6 4 v> R 

njj , #f/. 

/ — 8 (3 ■ S 0 6 sJo / 1 n 

c<o(« ) 

(i, k < ) ( ° • s') ^ 

4, _ — — — - 3 • i 2 5 *®/in 

2 (■•) ( S*) 

(a.,) Spring 2 ins.de spring 1; 1c , + 1< *- = SS ' 95 

(l) spring 2 on top 0 / spring 1: J- = ^ 4- 



tt/r 



^ 4c,+ 1<z S6.8o^ 4 3-125 



n. 
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x 2 = ^ Sin 3 0 

X| " y CoS 3 0 



t-e.. 



% 1 4 *<+ 



+ i * 2 y 1 s'" 1 



i ** 



4 



J OS' e-^Xic^.o 6 } = 2 . z , 7295 

1 2* loo 

-(<2= dlii = ? (7 1 - 6-5 1 ) (30 *'/) ro « #,/ 

2 a ^ -L_ L == 2-12 OSS X lO v»f 

*2. 75 

•*. = _L (2.237235 *,/ ) + '(2'|2fl?S X l o* ) 

u 4 

= 2* 2S 3 116 2 5 *«/ 'V.Vi 

SlYn'ilcuriy^ irke £^u ivalen t Jcumpm^ ConS-fcanf can te found axJ 

(uSiaj e.^ui va/Unoe of meti'c energies): 

| C.+ i- (••<.) = 0.37S 



! 4 7 ,j S'ta.mlesy sieei: e =■ 30 xi/ Ib/m 7 ' > G= n-sr x to &>/in 

For eack tute; 

X> = o • '•b o ;/ , J = o- 29", / = so' 




// 



Axfa,l stt'ff ne ss = A E _ tt / ^ ^ \ E 

J2 ” 4 ^ i 

= TL (»'30 1 -o.2g a ) p a 5 .^ l< ’ 6 ^ =2-730' 316 = * 
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Tors.oncU st;£fn ess = HJL 

32 l 





_ -n-Cii-Sxio 6 ) o 


•23*) x 2 




32 (Jo) 




For 


tie cut: e*ck center with 


6 tubes: 


Ax i‘<xL stij^£riess = £ - 1 


6,, 6 8 l * 8H(o 


TorS\ or»<a^t sti‘f,fr>e?S — g 


- 139-165 




f 1.48i] 



Assume sma.lt a_r\c,tes 9, «.nd Q z ) 9 Z = ^ tfj 9 { 

X - ho r« zoni&J- d\'£frla-ce.rr\e.Yyb of €•(** of ^i| - &\ 

1 " d * U Q f / L, 

vertrc cut cL's'pl'Z'Cemeni of C'Gj. of rncuts m 2 - r 2. - '* 1 z / 

= horijonta.1 drs pUeermen-fc of fc, cund h z - (r, + J?,) 

3 Z = Ver -fci’ca.1 <J.*s-f>kcemen£ Sf>rCn^t fc 3 ^ = &z i z = f>,X t 0,/f> z 



Eflui'v’&leYice 


of •< 


< .'metre emerges ^iVeS 






i 


- J. 
“ z 


^ ( 0 > )*+i J ife) Z+ 2 : 


m,(x, N ) 2 ' 


+ Z W 2 ( X 2 ) 


**• Tag - 


cT, + 


^z (V h.) 2 + n 2- 


+■ m z ^ 


( h/bxt 


EguiVculence 


o/ 


potent i'cU enerji'es <jj 


iVer 




i S = 


: i * 


/z ^ + t ^ \ + i 


*£, ©i 2 " -i- ■£ 0 2 - 




With 


^ 1 2 ~ fcf ’f’ » 


^34 " 


k 3 * 4 / ^3 + 



- &l (^l + ^ l) > ^2.~ /’j ^ ^ / b Z ^nel ©z ~ ^ A ' 

feej = (*' + **) (*' + *') + (■!$*;) ^ + ktl+ 7? ’ 





1.49 



9= f 



x, - *■ ^ 

1 b 



From eguiVaJence of fan etfc energies 

2. >i • 2. << • 2. 



t w «g = 7 



+ 4 rn_ x 






(-^) Z + (t) 
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1.50 



Let 0. s ounyula-r velocity of the motor ( input') 
Angula-r velocities of different jea.r sets cure : 



J moicr > ^1 J ^Zj ^3 
1 



, r r 



Equivalence of kinetic energies jives 

i ^ Q. = 1 cT cier e . + ± ^ + i *i ca J &l ea _j 

'* ^ - ( J rtoior + ' 7 f)+ ( T t+ + J %)( , ‘nz 

+ • - ' •+■ C ^2 N 






X m r>3 

m n^" 



J* 

2/^ 




1 rtC* ' 

1 . 


/ vu 
l Ol 


n 2M-1 \ 
n iNy 


IcclJ 


. z 

& U a .Jl 


- 


HI 


n 3 Y 2 


5 A 


n 2 


J 


n 2N -1 


\ 2 


n 


2N 


) 




1.51 



feuiVeulcnce of Kinetic energies yves 

I z H el - i 3-1 el + i wUre ^ = 

*t“ 31 + * 




1.521) 



When point A moves by distance x = x t , the walking beam rotates by the angle 
- -j-- 

3 • • a P - Xh ** 

This corresponds to a linear motion of point B: x B — y b ^2 7 

3 

and the angular rotation of crank can be found from the relation: 



X B 



= r c sin 9 C + cos 4> = r c sin 9 C + ^4 l - sin 2 9 C 

For large values of €4 compared to r P and for small values 

of x and 9 C , we have 

. a a X B x b A 

x B ~ r c sm 6 C = r c t) c or V c = — p 

r c c 3 r c 

The kinetic energy of the system can be expressed as 

T = ^ m h xj + J b 9 h + y J c K 

1 .21 .2 

Equating this to T = — nieq x — — x^, we obtain 

* p ^2 2 
€ 2 



, Jb 

m^q = m h + -2- 
^ 3 



+ Jc 



"3 r c 




'B 
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X 

@ Wlien mass m is displaced by x, the bell crank lever rotates by the angle 6 h = — . This 

makes the center of the sphere displace by x 3 = 6 h <? 2 . Since the sphere rotates with out 
slip, it rotates by an angle 



0 , = — = 



0 b it 



A r. 



1 r s 



The kinetic energy of the system can be expressed as 



T 



- z 






since for a sphere, J s = f m, * Equating this to T = ± m, i*. »e obtain 

, , 5 n 






“ m + J ° 7f + 5 1113 e\ 
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p c * = AcLmpinj force of Ci = ( *2. ^ 1 ) J 1 

f - da-mping force of c eg = Ce £ ( x 2 ” 
^ = F,+ F 2 .+ F 3 



(*>) 



= c, + c 2 + C 3 



K, ^ 






- 

Pee 

Ce & 

Si'nce 



*-/ = Xt 



■e% 

X 3 f kj — *2 + X; 



Fi - *0 

F 2 = <=2 " *2) 

Fj = C 3 



<•‘ 2 . 



F< 

c ( 



Peg - ~ F 2 = F 3 ' 



I L 4- — + ~ 

— = T, C 2- C 3 



i i nee >"eg - » i - 1 2 - ' 3 ' c eg, 

(c) eneryrej m a, 

7T CJ xf ^ Trc,09 A , Z + TT c i X* + 7T C 3 69 *3 

J _ i \ * ^ Q rk m J X? 



C c „ - 






l - 

wliere X , - e / 4 , X 2 = 0^2. a ' ric ^ X 3 - >9 3 

“ c ' + c * (rf* 

. - i-L- energies Jl i is ' p <*-t e.a m a, c-yole , 



TT C, 



Cfe ^' 



z 

&, = TT 


c ti 


z 

CJ &, -f T 


C t 2 65 0 2 + 


TT 


C -63 ^ 












/ ni \ 


l*/kere 


© 2 - 




Our\c\ 6? — 


0 I 


i VnJ/ 


= C -tl + 


c tz 


+ 
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© Damping constant desired = c = 1 lb-sec/in, viscosity of the fluid = 

H = 4 n reyn = 4 (10~ S ) lb— sec/in 2 . 









C = M ■ 


2d 1 

3 ir d s e (i + 




CO 

id 


* D * 


l 1 Assuming x = D/d as the unknown with € — 2 in, 



Eq. (1) can be written as 



c = M |3^lj(l + A) or 1 = (4 (10- 8 )) (-^-) x 3 (1 + | ) 



( 1 ) 



( 2 ) 



This gives x 3 + 2 x 2 — 53,051.52 — 0 , 

Using a trial and error procedure, the solution of this cubic equation can be found as 

36.92. Using D = 3 in, we get d = 3/36.92 = 0.08126 in. 
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1 . 58 ' 



-/*{ 



3 7T D 3 JL 



(l+i i)}; 



4 J 3 

yur: 45 yu reynolds 

Lerfc 4 = 0- 00\" > D = 2*4* 



an 4 



I0 5 - (4^,o- g )f ^ L (, 

t 4 (o.ool) 

1 - O’GS 17 " 



J — Ai'c^m e.ter of jpiston 
£ - clxCcJ Unj-kk of piston 
d = ra,cli'oJ cUa.ra.nCLe 



ajoove. e^u.a-tion ^iVef 
2. x o • O O I 



2« *+ 



1-41 



PRELIMINARY VERSION 12/14/2010 



© 201 1 Pearson Education, Inc., Upper Saddle River, NJ. All rights reserved. This publication is protected by Copyright and written permission 
should be obtained from the publisher prior to any prohibited reproduction, storage in a retrieval system, or transmission in any form or by any 
means, electronic, mechanical, photocopying, recording, or likewise. For information regarding permission(s), write to: Rights and Permissions 
Department, Pearson Education, Inc., Upper Saddle River, NJ 07458. 




1.59* 



To.nyeniiod velocity of >'»ner cylinder = -J & 
For sma-LL d, rode of cUnye of velocity of 

fluid is _ V.C 3 

dr ” "tr 

skea.r siress beiu/een cylinders is 

TT 

cund sh ea_r force is 




~K ^ D ^ ca ( i— A) 
F = r • Are^ - tr 7TP (l- A) = £“d 

Torque developed = Al t1 = F • 

For frncolL h > route of chcLmye of 

Velocity of fluid in ver-tfcod direcii'on. is 

dv _ r oa 

Ay i* H- res 

Siecur sir ess is V - f* dy ~ A 
Force on or ecu cl A = dF = dA 



7 orgwe hebuJeen boHom surfo-ces of cylinders is 

Mt2. = ff w/,ere 

„ r D ^ / , iTr yUdPTTP 

J / r 3 .dr d© = ^4 A 

r = o e = 0 , * 

TTmD 3 C 3 (l-L) . V W]> 

To-fca.1 -torque = /Y) t = Mti + - — L T + 

4 - a ^ 

Ex j>re?Wnj M.J- cuS C -( . v ■= 05 , we douxn^iny consicunt: 




r = 



/*JL-g dr dG 

U 



**« ^ 



c*-*) 



Z d 



+ 



TT 



A B ' 



32 A. 
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1.69) 



F = a x + b x 2 = 5 x + 0.2 x 2 
dF 

F(x) ~ F(xo) + — I io (* “ *o) 
dx 



At Xo = 5 m/s, F(xo) = 5 (5) + 0.2 (25) = 30 N , ^ | i„ = (5 + 0.4 i) | , - 7 and hence 

dx 

F(x) = 30+7 (x — 5) = 7 x — 5. . . . _ _ # 

Thus the linear ized damping constant is given by F(x) ~ 7 x — Cgg x or Cgg — 7 N— s/m. 

D am ping constant due to skin friction drag is: 

c = 100 [X P d 

Damping constant of a plate-type damper is: 

/i A 




1.70) 



( 1 ) 



c > = h 



( 2 ) 



where A = area of plates and h = distance between the plates. If the area of plates (A) 
in Fig. 1.42 is taken to be same as the area of the plate shown in Fig , 1 .107, we have A 
£ d. Equating (1) and (2) gives 

100 n € 2 d = ( 3 ) 

h 

from which the clearance between the plates can be determined as h — ^qq g '• 



© 



c = 



6 7r n € 






2 2 
a — r 



h 

a 

2 



-h 



When u = 0.3445 Pa-s, £ = 0.1 m, h = 0.001 m, a = 0..02 m, and r = 0.005 m: 



c = 



6 7T (0.3445) (0.1) 



(10 -3 ) 3 



(0.02 - 0.0005) 2 - 0.005 2 
= 4,205.6394 N-s/m 



0.02 2 - 0.005 2 
0.02 - 0.0005 



- 0.001 
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Ba,?(C deltas : f s to cm, - 0 > \ cm, a. = z > f = < 3 . 5 cm , 

yu - o- 3445 

CO n S"£ CLY\"t U/i’tA i ^CLSt'c dcctaj 

C - 4j 2 °5*623C> N-S/'m 

(Ou) r chcunjed to 1 cm ; new c= 2,617.7920 N - */m 

(t) ctiCLVi^eJ to Cm; new c- 3?, <5 £0 • 81 I 0 N-S/m 

(c) CL cka-n^ed £0 4 cm; new c = 3 8 , 7 5 4 • 8 60 N - s/m 
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1.75) 



* = 5 +■ 2 < 

A cos © = 5 
A Sin 9 - 2 



19 




1.76) 



- A e"' - A cos 9 -f *- A tin 9 

h = J(A CoS 9) l + (A An e) 2- = /?*+ 2 1 = ?*3852 

, "?2. = 3-4 i = b± + b z L 



Xj = If 2 t + <X 2 ^ , 

X = 4 . X 2 = ( Ou i _+ tl) -+- c (O'Z’h ^z') - 4- - 2 i 

- A = A cose + i A sin © 

A = <!*’■+ (-O* = 4-4721 
0 = = - 26 - 5651 “ 




1-77 



Zi = (3 - 4i), Z 2 = (l + 2 i) 
z = Zj — z 2 = (3 - 4i) - (1 + 2i) = 2 - 6i = A e* 6 

where A = 2 2 + ( — 6) 2 = 6.3246 and 6 — tan 



-l 



-6 



= tan 1 (—3) = — 1.2490 rad 




1.78) 



Zi = 1 + 2i, z 2 = 3 — 4i 

z = Zl z 2 = (1 + 2i)(3 - 4i) = 11 + 2i = A e 1 9 

where A = Vll 2 + 2 2 = 11.1803 and 6 = tan -1 (2/ll) = 0.1798 rad 




1.79 z = 



= 1 + ? 1 = C 1 + 2 W + 4 *3 . = ~ 5 + 10 1 = -0,2 + 0 .4i = Ae i 



z 2 3 — 4 i (3—4 i)( 3 + 4 i) 
where A = V(;-0.2) 2 x + (0.4 ) 2 = 0.4472 
-0.4 



25 



and 6 = tan 



-l 



0.2 



= tan -1 (-2) = - 1.1071 rad 
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1.80 



x(t) = X cos oj t , y(t) = Y cos (u> t + p) 
x 2 y 2 

(a) — = cos 2 0Jt, = cos 2 («t + p), 



X y 
XY 



COS (j> = 2 COS OJ t cos (oj t + p) cos p 

= cos 2 OJ t + cos 2 (ojt + p) —2 COS OJ t COS p cos (oj t + p) (1) 

Noting that cos 2 a = (1 + cos 2 a), Eq. (1) can be rewritten as 

4 + 4_ 2 2EL C os^ 

X 2 Y 2 XY ^ 

= — + — cos 2 oj t + -T + T cos (2 oj t + 2 p) — 2 cos oj t cos p cos (a; t + p) 
z z z z 



= 1+ ? 



_ 2ojt + 2ut + 2(b 2 oj t — 2 oj t — 2 d> 

2 cos cos — 



— 2 COS OJ t COS (f) cos (oj t + <f>) 

= 1 + cos (2 OJ t + <f>) cos <p — 2 cos OJ t cos (p cos (oj t + (p) 

= 1 + cos (2 oj t + <p) cos <p — 2 cos p j-T |cos (oj t + <p — oj t) + cos (a; t + p + oj t)|| 



= 1 + cos p cos (2 oj t + p) — cos p j cos p + cos (2 oj t + p) 

= 1 — cos 2 p = sin 2 p 
(b) When p = 0, Eq. (2) reduces to 

xy 



( 2 ) 



2 2 
x 2 Y 2 XY 

X 



x 

T 



y_ 

Y 



= o 



which gives X = ±— y. This indicates that the locus of the resultant motion is a 

7T 

straight line. When p = — , Eq. (2) reduces to 

2 

X 2 Y 2 

which denotes an ellipse with its major and minor axes along x and y directions, 
respectively. When p = 7r, Eq. (2) reduces to that of a straight line as in the case 
of p = 0. 
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Equation for resultant motion: 

I*- + - 2 yy cos2 * = sin2 * M 

x 2 

When y = 0, Eq. (l) reduces to — = sin 2 4> and hence: 

x = ± X sin <f> = ± 6.2 = OS in figure ( 2 ) 



|i=sii 

Y 2 

y = ± Y sin 4> = ± 6.0 = OT in figure 



When x = 0, Eq. (1) reduces to = sin 2 <f> and hence: 



It can be seen that 



OR = X cos <p = 7.6 in figure 



(3) 

(4) 




OS _ X sin> =tajl(t> = — = 0-8158 or 4> = 39.2072° 

OR X cos <j> 7.6 

From Eqs. (2) and (4), we find 

X = V(X sin + (X cos 4 >? = V(6.2) 2 + (7.6) 2 ” = 9.8082 mm 



( 5 ) 



Equations (3) and (5) give 



Y = 



6.0 



6.0 



sin (p sin 39.2072 c 



= 9.4918 mm 
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1.82 



(a>) x (t) = _A — cos f so t + <* \ 

mnfl ' ' 



1000 

*(*0 = i^o cos ^ = °’° o3 

■x(o) = - A jin - i 

<0 00 



where A i S in mm (E 

A COS o( - 3 (e z ) 

A Sirt o< = -20 (e 3 ) 

A = {(A coS 1 + (A sin o< )■ ^ ■= 20 • 2237 mm 

- tan“ 1 ( A ^ in } -ta.n _i (- 6* = - 81* +692°= - i-4219 rad 

\ A COS o( J \ y 

x(t) = 20*2237 COS (s-O-fc - 1.4-219} mm 

(k) COS(A-hB) = COS A coSB - Sin A Sin B 

Eg*(Ei) cam be expressed as x(t") = A cos sot-cos - A sin 5 o£.sm<< 

= A\ COS cot 4 - A z s *’ 0 09 -fc 
where 09= so, A 1 = A cos^, A^-Arin^ 
x(t) = (3 CoS 5ot +20 sin So t ^ mm 




1.83) 




1.84) 



x(t) = At COS cst ■+ A Z Sir. t 

d£(t ) — _ ^1 09 sin 03 1 4- A 2 09 COS cat , = - A COS <ot - A Z 09 2 Sin 03 1 

tit "t 

d- 2- * _ _ ^2 where os 2, is a, constant 

Hence x£t) is co Simple harmonic motion. 

(cu) Using trigonometric relations : 

X A (t) = y (cos 3 1 cos 1 - sin 3 1 Sin 1 } 

x 2 (t) = 10 C coS coS 2- _ s,n Sin 2 ) 

X (t) = x 4 (t) + = COS 3 1 (s COS 1+10 cos 2 )- sin 3t (5 sin 1 + 10 sin 2 ) 

If x(t) = A cos (cst 4 - <*) = A CoS csi: COS o( - A sin cot Sin °C , 

09 = 3 , A CoS <X = S coS 1+10 cos 2 = — 1* 4-599 , 

A sin e< = 5 *in 1 + IO Sin 2 = 13* 300 3 

A = y C A COS or ) 2 4 . (A Sin OC) 1 ’ = 13*3 802 

* _ tar) -i ^ A Sin^oc ^ _ tan -1 (-9- H04) = 96*2640°= 1*68 rad 
An 3 le between x t (t) and *(t) * 96* 2 6 4 0° - 5 7. 3 ° = 38 * 9 64° 

(t>) Us ing vector addition ; 

Por an arbitrary value of 
(c 9 t+i), harmonic motions 
Xi(t) and x^(t) can be shown 

as in the figure. From 

vector addition, we find 

x(t) d 13*38 CoS (c 9t + !• 68 ^ 
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(c) Using complex numbers: 




1.85) 



t nu.mcc.1 . i . \ , 

r i (Mt +!') ) D P f 5 / ( ^ M 
*!« = Be {A e \ = Re i S e 



i(u 5 -fc + 2 .) 



} 



x 2 (*) = Re {At + e 

if x(t) = Re { A e i <' 6st+ • 

A coS(3t + «i) = Ai cos(3t+l) + A 2 cos (it + O 
r.e. A (cos it CCS * - Si-' 3t sin <)=5 (cos 3t- COS 1 - sin it- Sin l) 

+ 10 (cos 3t- CoS 2. — Sin 3t • Sin l') 

re- A COS ^ = 5 cos 1 + 10 CoS z , A sin <V = 5 Sin l+IO Sin 2 

A= 13-3 802 , * = 1-6* raA , 

r * (3t + 1*^8) ^ 

x("t) = Re £ 13 *3 802 . 6 J 

*(0= 10 sin (os-t + GO°) - *i(t) + x z(t) 

where X i (t') = 5SiV(«t+30 5 ) a-nd X z (t> = A S.n (cst + << j 

,0 (sin cat COS 60“ 3- COS cat sin «0* ) = 3 (Ain cat Cos 30* + Cos cot sin 30*) 

+ A (Sin CSt-CoS =*°+- COS Csi-S'tn Of ) 

(0 cos <oO° = ? CoS 30° + A CoS o<° ; A cos A 0 = 0* 6699 

10 sin 60° = 5 Sin 30°+ A sin of 0 ) A Sin = 6-/603 

A = 7°‘ 6699 2 + 6-I6 0 3 2- = 6.|?G6 



c< = to.n _1 ( 6,,6o3 /o. G6 99 s ) = 83-7938 

= 6*19 66 sin ( uyt -j- 8 3*7938 ) 




.*. is not harmonic 
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1.88 



x(t) = i cos - cos nt 

- Ilf cos JL t + 7t 2 ‘ cos 7rt 

s ^ 

i’S not harmonic 



— consla-rvfc himes *00 




(1.89 



x(t) = *i(t) ■+■ = 3 s,n 3ot + 3 s ' n z9t 



Sin A 4- Sin B = 2 Sin 



A + S 



cos 



A-B 







(1.90) 



Since 

x(0 = (g cos 

This e<fruaJ:ior) shows 
hhod -the amfjLMruele 
(G cosi|.) varies with 
time between co maximum 

Value of 6 and a* 
minimum value of O • 

The frequency 0 f this 
oscillation ( beat 
f re^uency^ lS 03 y- 1. 

Note; Becct frequency is twice the freguency of the terrn^ 

6 cos t since two pea-KS bass i n each cycle of (<s 

2 — 

The resultant motion of two harmonic motions having identical amplitudes (X) but 
slightly different frequencies (co and co + <5co) is given by Eq. (1.67): 




rnmimum 

am p litudes 



x(t) = 2 X cos 



COt + 



5co t 



cos 



5 co t 



Thus the maximum amplitude of the resultant motion is equal to 2X and the beat 
frequency is equal to <5co. From Fig. 1 . 113 , we find that 2X ~ 5 mm or X = 2.5 mm 
and 

&0 _ 2 7T _ 2 7T _ 2 7T 

T ~ r beat ~ r larger 2 (12.6 - 4.2) 



= 0.374 rad/sec 



or 60 J = 0.748 rad/sec and co + 



S co 



2 7 r 



2 7 r 



= 6.2832 rad/sec 



2 ^smaller t 

Hence co = 6.2832 - 0.3740 = 5.9092 rad/sec. Thus the amplitudes of the two 
motions = X = 2.5 mm and their frequencies are co = 5.9092 rad/sec and co + <5co 
= 5.9092 + 0.7480 = 6.6572 rad/sec. 




1.91) A = 0-05 m 



C3 — IO 



Hg- =. 6 2*832 ra^d/sec 



period = T — 3JL - LZI — _ 0 

r 05 6 2*832 



sec 



■maximum velocity = A 05 = 0*05 x 62*832 — 3-14(6 rn / s 

Uera,h'or> = A OS 2 = o*o5 (62*832) 2 = I <3 7*393 Vs 2 



maximum a^ccet 
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(1.92 




^ - is c(>s = 94-248 ra A/ sec 

~ = o.S* = 0.5(9-81) = 4-905 m/s* = Ac* 

A= amplitude = 4 - 905 /( 34 . 242 / - o-ooo552Z - 

= ma-x. vetoed = A 09 = 0-05204 ™/± 

JC - A COS oat , - A = 0-25 mm , X = - A CoS ^ ^ 

A69 i = 0 . 4 ^ = 392 A J 09 2 = 2s2 -Va = 1569 6 (ro- /s) 

Ob£ro ^ o£ bun** .25-28.7 roJ/s = 
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2 7T t , 

x(t) = X sin — — > 



x rms 



t J* 



•2 2 ^ 7T t ^ + 

1 sun cLt 



1 — cos 



T 

4 7T t 





For even functions, = Vz 

in - * "*»* M 



= A f" f *(£') s, ' r * ncat-Jfc + J sin ncsi-dt 

L - X ° (Ei) 
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with x(t) denoting an approximation to the exact x(t) given by Eq. (1.70). The 
error to be minimized is given by 

( 2 ) 
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E = / e 2 (t) dt 
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where e(t) = x(t) — x(t) (^) 
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Rearranging Eq. (4) gives 
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0 for mrD 

tt/u 

J x(t) cos n <jJ t dt = 

— 7t/u 



( 6 ) 



a* 7T 



for m = n 
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In a similar manner, we can derive: 
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It can be observed that Eqs. (9) and (10) are simlar to those of Eqs. (E.3) and 
(E.4). 
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23 'a C7 c ' -22.00 i -21.25 5.69 ',-19.05 11.00 ■ -15.56 15.56 

24 ! o-60o ' 0.00 [ 0.00 0.00 j 0.00 0, _ 00 j_ 0 -°_°__ 0 J° 0 - 

239.00 -241.90 282.30 39.72 147.18 45.26 -4.88 

19.92 -20.16 23.53 3.31 12.26 3.77 -0.41 



% 

%Programl .m 

%Program for calling the subroutine FORIER 



) 



24 

** * 2 .* 

kP< 
- 1 



% 

%Run "Programl .m" in MATLAB Command Window. Programl.m and forier.m should be 
%in the same file folder, and set the path to this folder 
%Following 6 lines contain problem- dependent data 
n=16; 



m=3 ; 

time=0 . 32 ; 

x=[9 13 17 29 43 59 63 57 49 35 35 41 47 41 13 7] ; 
t=0. 02:0. 02:0. 32; 

%end of problem-dependent data 
%Following line calls subroutine forier.m 
[azero, a,b, xsin,xcos] =forier (n,m, tim e,x, t) ; 

%following outputs data 

fprintf ( ' Fourier series expansion of the function yL(t)\n\n'); 
f print f ( ' Data : \n\n ' ) ; 

fprintf ( 'Number of data points in one cycle = %3.0f \n',n); 
fprintf ( ' \n'); 

fprintf ( 'Number of Fourier Coefficients required = %3.0f \n',m); 
fprintf ( ' \n'); 

fprintf ( 'Time period = %8.6e \n\n' 7 time); 
fprintf ( ' Station i ' ) 

fprintf ( 'Time at station i: t(i) ; ) 

fprintf ( 'x(i) at t(i)') 
for i=l:n * 

fprintf ('\n %8d%25.6e%27.6e ' , i , t (i) # x (i) ) ; 

end 

fprintf ( ' \n\n' ) ; 

fprintf ( 'Results of Fourier analysis : \n\n' ) ; 
fprintf ( ' azero=%8 . 6e \n\n' , azero) ; 

fprintf ( 'values of i a(i) b(i)\n'); 

for i=l:m 

fprintf ('%10.0g %8.6e%20.6e \n' , i, a (i) ,b (i) ) ; 

end 
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% 

%Subroutine forier.m 



function [azero , a, b, xsin, xcos] =f orier (n, m, time, x, t) 
pi=3 .1416; 
sumz=0 . 0 ; 
for i=l:n 

sumz=sumz+x(i) ; 

end 

azero=2 . 0*sumz/n; 
for ii=l : m 
sums=0 . 0 ; 
sumc=0 . 0 ; 
for i=l:n 

theta=2 . 0*pi*t (i) *ii/time; 
xcos (i) =x(i) *cos (theta) ; 
xsin(i) =x(i) *sin( theta) ; 
sums = sums +xs in (i) ; 
sumc=sumc+xcos (i) ; 

end 

a (ii) =2 . 0*sumc/n; 
b (ii) =2 . 0*sums/n; 
end 

» programl 

Fourier series expansion of the function x(t) 



Data: 



Number of data points in one cycle = 16 

Number of Fourier Coefficients required = 3 

Time period = 3 . 200000e-001 



l i 


Time at station i: t(i) 


x(i) at t(i) 


1 


2 . 000000e-002 


9 . 000000e+000 


2 


4 . 000000e-002 


1 . 300000e+001 


3 


6 . 000000e-002 


1.700000e+001 


4 


8 . 000000e-002 


2 . 900000e+001 


5 


1 . 000000e-001 


4 . 300000e+001 


6 


1.200000e-001 


5 . 900000e+001 


7 


1 . 400000e-001 


6 . 300000e+001 


8 


1 . 600000e-Q01 


5 . 700000e+001 


9 


1 . 800000e-001 


4 . 900000e+001 


10 


2 . 000000e-001 


3 . 500000e+001 


11 


2.200000e-001 


3 . 500000e+001 


12 


2 . 400000e-001 


4 . 100000e+001 


13 


2 . 600000e-001 


4 . 700000e+001 


14 


2 . 800000e-001 


4 . lOOOOOe+OOl 


15 


3 ..000000e-001 


1 . 3 00000e+001 


16 


3.200000e-001 


7 . 000000e+000 



Results of Fourier analysis : 



azero=6 ..9750Q0e+001 



values of i 
1 
2 
3 



a(i) 

-2.084870e+001 
-1 . 456887e-r000 
-5 . 402900e+000 



b(i) 

-3 . 915985e+000 
-1 . 144979e+001 
3.353473e+-000 
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% Exl_119.m 
for i = 1: 101 

t (i) = 0.32* (i-1) /100; 

x(i) = 34.875 - 20 . 8487*cos (19 . 635*t (i) ) - 3 .9160*sin(19 . 635*t 

- 1 .4569*cos (39 . 27*t (i) ) - 11 . 4498*sin(39 . 27*t (i) ) . . . 

- 5 . 4029*cos (58.905*t(i) ) + 3 . 3535*sin (58 . 905*t (i) ) ; 



i)) ... 



end 

plot (t,x) 
xlabel('t') ; 
ylabel ( 'x(t) ' ) ; 





% Exl_1.120.m 
u = 0.3445; 

1 = 10 ; 
hO = 0.1; 
aO = 2; 
rO = 0 . 5 ; 

% First case, r changes 
for i = 1:101 

r (i) = 0.5 + (i-1) *0.5/100; 

cl (i) = ( 6*pi*u*l/ (hO A 3 ) ) * ( (aO - h0/2) A 2 - r(i) A 2 )... 
* ( (aO A 2-r (i) A 2) / (a0-h0/2) - hO ); 

end 
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% Second case, h changes 
for i = 1:101 

h (i) = 0.05 + (i-l)*0. 05/100; 

c2(i) = ( 6*pi*u*l/(hU) A 3) ) * ( (aO - h(i)/2) A 2 - rO A 2 )... 

* ( (aO A 2-rO A 2) / (a0-h(i) /2) - h(i) ); 

end 

% Third case, a changes 
for i = 1:101 

a (i) = 2 + (i-1) *2/100; 

c3(i) = ( 6*pi*u*l/ (hO A 3) ) * ( (a(i) - hO/2) A 2 - rO A 2 )... 

* ( (a (i) A 2-rO A 2 ) / (a (i) -hO/2 ) - hO ); 

end 

subplot (311) ; 
plot (r, cl) ; 
xlabel ( 'r ' ) ; 
y label ( ' c (r) ' ) ; 
subplot (312 ) ; 
plot (h, c2) ; 
xlabel ( 'h' ) ; 
ylabel ( ' c (h) ' ) ; 
subplot (313 ) ; 
plot(a,c3); 
xlabel ( 'a' ) ; 
ylabel ( ' c (a) ' ) ; 
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% Exl_121.m 

for i = 1:101 

x(i) = (i-1) *4/100 ; 

ka(i) = 1000*x(i) - 100*x(i) A 2; 

kb(i) = 500 + 500 *x(i) A 2; 

end 

plot (x, ka) ; 
hold on 

plot (x, kb, ' — ' ) ; 
xlabel ( 'x' ) ; 

ylabel ( ' ka: solid line kb: dash line'); 





% Exl_122.m 

for i = 1:201 

t (i) = (i-1) *30/200; 
xl (i) = 3*sin (30*t (i) ) ; 
x2(i) = 3*sin (29*t (i) ) ; 
x ( i ) = xl ( i ) + x2 ( i ) ; 

end 

plot (t,x) ; 
xlabel ( ' t ' ) ; 
ylabel ( 'x' ) ; 
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Unbalanced force developed = P = 2 m <J r cos oj t, range of force 
range of frequency = 25 - 50 Hz = 157.08 - 314.16 rad/sec. 
Parameters to be determined: m, r, u. 

Let r = 0.1 m. To generate 100 N force at 25 Hz, set: 

= 100 = 2 m (157.08) 2 (0.1) 



= 0-100 N, 



max 



which gives 



100 



m = 



= 0.0202641 kg = 20.2641 g 



2 (157.08) 2 (0.1) 

To generate 100 N force at 50 Hz, set: 

P max = 100 = 2 m (314.16) 2 (0.1) 

which yields 



100 



m = 



2 (314.16) 2 (0.1) 



= 0.0050660 kg = 5.0660 g 



(l125) 




Goal: Weight to be maintained at 10 ± 0.1 lb/min . , 

Parameters to be determined: Angular velocity of crank (w), lengths of crank and 
connecting rod, dimensions of the wedge, dimensions of the orifice in the hopper, 
dimensions of the actuating rod, and dimensions of the lever arrangement. 

Given: Density of the material in the hopper. 

Select^ Abased on available motor. Determine the dimensions of the orifice in the 
hopper which delivers approximately 10 lb/min (assuming continuous flow of 
material ). For trial dimensions of the wedge, determine the increase/ decrease in 
the size (diameter) of the orifice. Choose the final dimensions of the wedge such 
that the material flow rate delivered by the orifice lies within the specified range. 

Force to be applied = 200 lb, frequency = 50 Hz = 314.16 rad/sec. 

Procedure : 



1. Select a motor that provides, either directly or through a gear system, the 
desired frequency. Assume that it is connected to the cam. 

2. determine the sizes and dimensions of the plate cam and the roller. 

3. Choose the dimensions of the follower. 

4. Select the weight as 200 lb. From the geometry, determine the range of 
displacement (vertical motion) of the weight. 

5. Determine the force exerted due to the falling weight. 
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Considerations to be taken in the design of vibratory bowl fee/ders: 

1. Suitable design of the electromagnet and its coil. 

2. Radius of the bowl and the pitch of the spiral (helical) delivery track. 

3. Tooling to be fixed along the spiral track to reject the defective or out-of- 
tolerance or incorrectly oriented parts. 

4. Design of elastic supports. 

5. Size and location of the outlet. 




Axial spring constant of each tube — k ^ • 

Let diameter of each tube be 0.01 m (1 cm) with thickness 0.001 m (1 mm). Then 



A = ZL (D 2 - d 2 ) = - (0.01 2 - 0.008 2 ) = 28.27 (10" 6 ) m 2 
4 4 



This gives 

k _ (28.27 (10 6 )) (2.07 (10 11 )) _ 29 26 (10 5 ) N/m 
2 

Since 76 tubes are in parallel, we have the total axial stiffness as: 
keq = 76 k = (76) (29.26 (10 5 )) = 222.38 (10 6 ) N/m 
The polar area moment of inertia of each tube is 

j = JL (D 4 — d 4 ) = ^ (0.01 4 - 0.008 4 ) = 580 (10“ 8 ) m 4 

Torsional stiffness of each tube is given by 

G J _ (79.6154 (10 9 )) (580 (10 ))^ _ 231 ^3^ N_ m / ra( i 

£ 2 

For 76 tubes in parallel, equivalent torsional stiffness will be: 
k t = (76) (231 (10 3 )) = 17.56 (10 6 ) N-m/rad 

“eq 
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